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ABSTRACT 

We report unique EVLA observations of SN2011fe representing the most sensitive radio study of a 
Type la supernova to date. Our data place direct constraints on the density of the surrounding medium 
at radii ~ 10 15 — 10 16 cm, implying an upper limit on the mass loss rate from the progenitor system of 
M < 6 X 10~ 10 M-o yr _1 (assuming a wind speed of 100 km s _1 ), or expansion into a uniform medium 
with density ncsu J$ 6 cm~ 3 . Drawing from the observed properties of non-conservative mass transfer 
among accreting white dwarfs, we use these limits on the density of the immediate environs to exclude 
a phase space of possible progenitors systems for SN 2011fe. We rule out a symbiotic progenitor system 
and also a system characterized by high accretion rate onto the white dwarf that is expected to give 
rise to optically-thick accretion winds. Assuming that a small fraction, 1%, of the mass accreted is 
lost from the progenitor system, we also eliminate much of the potential progenitor parameter space 
for white dwarfs hosting recurrent novae or undergoing stable nuclear burning. Therefore, we rule 
out much of the parameter space associated with popular single degenerate progenitor models for 
SN2011fe, leaving a limited phase space largely inhabited by some double degenerate systems, as well 
as exotic single degenerates with a sufficient time delay between mass accretion and SN explosion. 
Subject headings: supernovae: individual (SN2011fe) — supernovae: general — novae, cataclysmic 
variables — binaries: general — circumstellar matter 



1. INTRODUCTION 

Type la supernovae (SNe la) are luminous stellar ex- 
plosions that display remarkable homogeneity in their op- 
tical prop erties. C omprising ~30% of all SNe in the local 
universe (|Li et al.l 1201 laft . SNe la are exploited as bea- 
cons for cosmography and represent the dominant source 
of iron-peak elements; however, we still do not under- 
stand the nature of their progenitor system(s). There is 
a general consensus that SNe la mark the fatal disrup- 
tion of white dwarfs (WDs) near the Chandrasekhar limit 
(Mch = 1-4 Mq), and accretion from a binary compan- 
ion is required to reach sufficient WD mass and trigger 
the SN explosion (jHillebrandt fc Niemeverl 12000ft . How- 
ever, the nature of this binary companion remains un- 
clear, with potential progenitor systems falling into two 
broad classes: single degenerate (SD; containing a WD 
and a main- sequ ence, sub-giant. He s t ar, or red-gian t 
companion star; iWhelan fc Ibenl 119731: iNomotol Il982l ). 
and double degenerate (DP; in which two WDs merge; 
!Webbink|[l98l llben fc Tutukovi [l98l . In both scenar- 
ios, the evolution of the progenitor system shapes the 
circumbinary environment. 

Radio observations provide a sensitive probe of cir- 
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cumstellar material (CSM) surrounding SNe; when the 
Shockwave plows into this material, it accelerates par- 
ticles and amplifies the magnetic field, prod ucing syn- 
chrotr on emission that peaks in the cm band ( Chevali er 
1982b). The early synchrotron signal traces the CSM 
particle density, ncsM, on a radial scale of < 1 pc, 
the region shaped by the final stage of progenitor evo- 
lution. Radio emission is routinely detected from near by 
(d L < 10 Mpc) core-collapse SNe (jWeiler et al.l 12002ft — 
explosi ons that mar k the death of massive stars, M > 
8 Mq ()Smarttll2009D . However, radio emission has never 
been detected from a young SN la, despite searches 
spanning the past three decades (jPanagia et al.l 120061 : 
IHancock et ani20Tl . 

On 2 011 August 24 U T, the Palomar Transient Factory 
(PTF; lLaw et all I2009D discovered an optical transient 
withi n the nearby spiral galaxy M101 at a distance of 6.4 
Mpc (jShappee fc StaneHl20Tll Figure [J). Based on pre- 
and post-discovery imaging, the e xplosion date is we ll 
constrained to 2011 Aug 23.69 UT (jNugent et al.ll2011ft . 
Early optical spectroscopy revealed the tra nsient to be a 
SN o f the hydrogen-poor Type la class (|Nugent et al.l 
120 111) , the nearest such event discovered in 25 years, 
dubbed SN2011fe (PTFllkly). Detailed optical studies 
indicate that SN2011fe is p hotometrically (Nug ent et al.l 
120111 : [Margut ti et al.l [2012ft and spectroscopically (Par- 
rent et al. 2012, in preparation) a normal Type la. 

Here we present deep Expanded Very Large Array 
(EVLA) radio observations of SN2011fe, obtained just 
two days after explosion and enabling a factor of ten im- 
provement in sensitivity over previous limits. We note 
that radio observations of SN 2011fe were recently re- 
ported by iHoresh et al.l (|2012ft ; however, our data are 
significantly deeper and span a longer timescale. We pro- 





EVLA 5.9 GHz 

Fig . 1. — Images of the environs of SN 2011fe. Th e three-color i mage of the host galaxy, M101, features 21 cm emission (green; Walter ct al. 
2008) tracing neutral hydrogen, Ha emission (blue; Hoopcs ct al. 2001) tracing star formation, and 3.6 /an emission (red; [Dale et al.H2009l ) 
tracing the stellar mass. The bottom left image shows detail of the explosion site of SN2011fe, while the bottom right image shows our 
combined EVLA image at 5.9 GHz centered on the SN position (black circle has a radius of 10 arcsec). 
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vide a theoretical analysis of the data that and a fresh 
perspective on the relation between radio luminosity and 
CSM density around SNe la. We discuss our observations 
and interpret the measured upper limits in Sections 2 and 
3, and then consider the implications for a wide range of 
SN la progenitor scenarios and previous observations in 
Sections 4, 5, 6, and 7. 

2. OBSERVATIONS 

We initiated radio observations of SN2011fe on 2011 
Aug 25. 8 UT (At ~ 2 da ys after explosion) with the 
EVLA (iPerlev et al.l 12011ft under our program AS1015 
(PI Soderberg: IChomiuk fc Soderberdl2llill) . No radio 
counterpart was detected at the optical SN position to a 
limit, F v < 19 /iJy (3er; Figure[2]). We continued observa- 
tions of SN2011fe through our Director's Discretionary 
Time program 11B-217 (PI Soderberg) over six epochs 
that span three weeks following the explosion, as de- 
scribed in Table 1. Our subsequent EVLA observations 
reveal similar upper limits (Figure [5]) ■ Observations were 
carried out in A configuration using 2 GHz of bandwidth 
and recording four polarization products; one baseband 
of 1 GHz width was centered at 5.0 GHz while the other 
was centered at 6.75 GHz. Each epoch had one hour du- 
ration, yielding 40 minutes on source, except 2011 Aug 
26, when we observed for two hours. Data were cali- 
brated using J1349+5341 and 3C286, and were reduced 
using standard routines in the Astronomical Image Pro- 
cessing System (AIPS). Each sideband was imaged in- 
dividually; the 6.75 GHz image was then smoothed to 
the resolution of the 5.0 GHz image, and the two noise- 
weighted images were averaged together. We also cre- 
ated a stacked image by concatenating all uv data for 
a given sideband and producing a final image from the 
combined data set. Again, the two sidebands were av- 
eraged together in the image plane, producing an image 
that reaches the thermal noise limit at 5.9 GHz. Our 
stacked EVLA data place a strong constraint on the ra- 
dio luminosity of L v < 3.2 x 10 23 erg s _1 Hz -1 on an 
averaged date of 9.1 days after explosion (Figure [TJ. 

The resolution of our EVLA images is 0.56" x 0.39". 
The flux densities quoted in Table 1 are m easured at the 
SN coordinates given by iLi et al.l (|2011bft . The nearest 
point displaying any significant flux in our stacked image 
is a 4.5<7 (10.2 pjy at maximum) po int-like peak located 
0.7" away from the ILi et al.l (|2011bft position (at J2000.0 
RA = 14 h 03 m 05.799 s , Dec = +54°16'24.75"). However, 
when we image the RR and LL correlations of the EVLA 
polarizations separately, we find that no source is visi- 
ble in the RR image (3 ± 3 pJy in RR, as opposed to 
16 ± 3 /zJy in LL), thus implying that it is it is likely an 
instrumental artifact, and unlikely to be a real source. 
We also investigated if the position of the radio peak 
could possibly be consistent with SN2011fe by checking 
the world coordinate systems of our EVLA image and the 
H ST ACS /F814W image (GO-9490, PI K. Kuntz) used 
by ILi et all (|2011bft for consistency. We cross-compared 
the positions of two compact H II regions located ~1.5' 
northeast of SN2011fe that are visible in both the HST 
and EVLA images. For this comparison, we also made 
use of a narrow-band Ha image from _ffST/ WFPC2 (GO- 
5210, PI J. Trauger), as these H II regions are much more 
clearly defined in F656N than in F814W. We find that 
the coordinate systems are not offset by more than 0.2", 



TABLE 1 

EVLA Observational Data for SN2011fe 



Date Flux Density Time Since Explosion 



(UT) (pjy) (Days) 



2011 Aug 25.77 7.8 ± 5.8 2.1 

2011 Aug 27.72 0.4 ±6.1 4.0 

2011 Aug 29.98 -8.4 ±6.9 6.3 

2011 Sept 2.68 0.7 ± 4.3 9.0 

2011 Sept 8.04 -2.7 ±6.6 14.4 

2011 Sept 12.93 0.6 ± 6.0 19.2 

Stacked: Aug 25.8-Sept 12.9 -0.6 ± 2.2 2.1-19.2 



and thus the 4.5c radio peak is positionally inconsistent 
with SN 201 lfe. 

3. DENSITY CONSTRAINTS FROM RADIO 
UPPER LIMITS 

We make predictions for the radio emission from 
SNe la by analogizing them to other hydrogen-poor 
SNe which are detected in the radi o, namely Type lb 
and Type Ic SNe (SNe Ibc; e.g.. iBerger et alt l200l 
iChevalier fc Franssoiall2006t iSoderberg et al.ll2010ft . SNe 
Ibc are generally understood to mark the gravitational 
core collapse of massive stars that have been stripped of 
their hydrogen envelopes prior to explosion. The analogy 
between SNe la and Ibc is motivated by several common- 
alities, including (i) their compact progenitors, leading 
to mildly-relativisti c Shockwaves a t the time of break- 
out from the star ([Colgate! 11975 iNakar fe Sari 12012ft . 
and (ii) a common polytrope of 7 = 4/3, appropriate 
for both massive (relativistic) WDs and the radiative en- 
velopes of Wolf-Rayet stars. It is estimated that for SNe 
arising from compact progenitors, the outermost density 
profile of the ejecta (hosting the fastest moving material, 
which gives rise to the radio emission) is characterized 
by a steep power-law, psN °c r -10 2 ; (M atzner fe McKeel 
1999). The coupling of mass to ejecta velocity is sim- 
ply a function of the bulk explosion energy (-E51, in 
units of 10 51 erg) and ejecta m ass (M e j; IBerger et al.l 
l2002t [Chevalie r fe Frans son 2006). These parameters are 
well constrained in the case of SNe la; E51 = 1 and 
M e i = M c h = 1.4 Mq are fiducial values (Mazzali ct al. 
|200W> . 

The interaction of the Shockwave with the immediate 
environment is then described by a self-similar solution, 
such that the evolution of the Shockwave radius (R s ) and 
its velocity (v s ) are determined by the properties of the 
environ ment: its particle density (ticsm) and its radial 
profile (|Chevalierl 11982a!) . This dynamical interaction 
accelerates particles from the circumstellar environment 
into a power-law energy distribution, N(E) = No E~ p , 
where N(E) is the density of relativistic electrons at a 
given energy E, above a minimum energy, E m . We fur- 
ther assume that constant fractions of the post-shock en- 
ergy density is shared between relativistic electrons (e e ) 
and amplified mag netic fields (eg) and a shock compres- 
sion factor, 77 — 4 (Chevalier & Fransson 2006). In this 
framework, 

B = v^tt e B Pcsm v 2 s , (1) 
N = (p-2) e e p CSM v 2 s E^\ (2) 
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Fig. 2. — Deep limits on the radio luminosity of SN2011fe. Our 5.9 GHz upper limits (black arrows; 3cr) from six EVLA epochs 
are compared with model light curves for three different progenitor wind mass loss rates: M/v w = 1.5 X 10 — 8 loo^kn-Ts- 1 3 X 
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0.033 and 0.066) and assumes e e 



0.01 — 0.1 (with the intermediate boundaries between shades 
151 



100 km s~ 

corresponding to eg =0.033 and 0.066) and assumes e e = 0.1, M cj = 1.4 Mq, and E K = 10 bi erg. We also plot our stacked limit as a 
white arrow, used to constrain a model for uniform-density CSM with a partic le density of < 6 cm~ 3 (for eg = 0.1; dashed black line). 
Finally, overplotted as grey arrows are additional data from Horcsh ct al. (2012), scaled to 5.9 GHz. 



and 



E = 



(p - 2)e e nm p v 2 s 
(p - l)r)(n e /ni) ' 



Here fj, is the mean molecular weight of the CSM as- 
suming it is neutral, and (n e /rii) is the number ratio 
of electrons to ions in the shocked gas (for solar abun- 
dance, [i = 1.4 and n e /rii = 1.14). The shock-accelerated 
electrons gyrate along the magnetic field lines, giving 
rise to synchrotron radiation that peaks in the radio 
shortly after explosion with a spectrum, L v oc here 
(3 = — (p — l)/2. The high velocity Shockwaves of Type 
I SNe minimize opacity to external free-free absorp- 
tion, such that synchrotron self-absorption dominates the 
spectral energy d istribution, supp ressing the emission at 
low frequencies (IChevalierlfToM ). For sources charac- 
terized by synchrotron self absorption, the shock radius 
may be robustly inferred from the synchrotron light curve 
(Rcadhcad 1994), thus enabling a direct tie between the 
observed radio p roperties and the SN self-similar solution 
(|Chevaliedll998l) . 
For SN2011fe, we assume p = 3, e e 



and a r ange of 6» = [0-01 — 0.1] as inferred for 
SN Ib/c (IChevalier k , Fra"nssoil[20Q6l : IZhang et al.ll200l 
( 3 ) ISoderberg et all 1201 ID . In this framework, the radio lu- 
minosity at a given frequency and time after explosion 
dep ends only on nc, 9,y\ , and is calculated using Equation 
1 of lChevalierl (|1998f ) Therefore, radio observations probe 
the environment of SNe la and may distinguish between 
different progenitor systems (see Chomiuk et al. 2012, in 
preparation, for a detailed discussion). 

3.1. Wind Density Profile for the CSM 

A commonality among SD progenitor channels is the 
accretion of material onto a WD due to its interaction 
with a donor star. While poss ible progenitor scenarios 
span a wide parameter space (fHillcbra ndt fc Niemeverl 
2000), most scenarios may be characterized by the WD's 
accretion rate, M acc , and efficiency of accretion, e acc , 
such that the rate of mass lost into the local e nvirons 
is M loss = (1 - e acc )Af acc (INomoto et a lJ l2007D . Here 
we define, ei oss = (1 — e acc ). In this non-conservative 
0.1, mass transfer scenario, we assume that the lost material 
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is blown outward with a wind velocity, v w , and radially 

distributed as rtcsM = M / (Airv w r 2 \ 

Using the similarity solutions of lChevalierl (|1982al ) and 
the SN la ejecta density profile of Chomiuk et al. (2012, 
in preparation), we estimate the Shockwave radius to be: 

R s = 5.5 x 10 14 (M /v w )~ - 12 £ 5 Y 3 M- h - 31 t% S8 cm, (4) 
where t\o is the time since explosion in units of 10 days. 
For t w = 0.21, M/v w = lO" 10 ^ P^ si = 1, and 
M c h = 1, we find R s = 2.2 x 10 15 cm and the Shockwave 
velocity (v s = 0.88R s /t) is 0.35c. 

When the SN Shockwave interacts with such a medium, 
the radio emission is suppressed at early times by syn- 
chrotron self absorption, and declines at later times due 
to the decelerating Shockwave and decreasing density 
of the CSM. With the radius evolution and assump- 
tions described above, we calculate a family o f radio 
ligh t curves using the fo r malis m of iChevalieil (|1998f ) 
and iChevalier fc Franssonl (|2006T ) . In the optically-thin 
regime (after light-curve peak, which is valid here; Figure 
[2]), for a given luminosity limit, our constraint on M /v w 
scales as e^ ' 7 . 

Our measured upper limits on L u measured At ~ 2.1 
days after explosion provide the deepest limits on the 
density of the environment surrounding SN 20 life to 

date, corresponding to M/v w < 6 x 10~ 10 1 ^ g> k ^ r s _ i 

(assuming e B = 0.1) or M/v w < 3 x 10~ 9 ™ \^ s -i 
(assuming es = 0.01) for the SN2011fe progenitor sys- 
tem. These l imits are a factor of ~ 20 deeper than those 
presented bv lHoresh et al.l (|2012D . due to our calculation 
from first principles of v s and our self-consistent treat- 
ment of E m . 

Mass loss from a SD progenitor may be concentrated 
toward the equatorial or polar re gions of the orbit (e.g., 
iMohamed fc P odsiadlowski [20111 ). Assuming that the 
asymmetric mass loss still follows an ncsM °c f 2 ra- 
dial profile, and covers a solid angle VL (units of ster- 

radians), the radio luminosity scales as L v oc (jM 
The increase in CSM density due to asymmetry wins out 
over the slightly stronger deceleration of the Shockwave, 
leading to an increase in the radio luminosity for a given 
M/v w , if material is expelled asymmetrically. Therefore, 
our assumption of spherical symmetry is conservative. 

3.2. Uniform Density Profile for the CSM 

We also consider a uniform density medium, as might 
be expected if the SN is exploding into the ambient inter- 
stellar mediu m. In this c ase, ag ain using the similarity 
solutions from IChevalier] (|1982a[) and the density profile 
of SN la ejecta from Chomiuk et al. (2012, in prepara- 
tion), the SN Shockwave radius evolves as: 

R — 1 A v ml6 „-0.10 ttO.35 jljt- 0.25 ,0.71 fr\ 

R s - 1.4 x 10 n CSM h 51 M ch t 1Q cm. (5) 

For t w = 0.21, n C SM = 1 cm- 3 ,£ 51 = 1, and M ch = 1, 
we calculate R s = 4.5 x 10 15 cm and a Shockwave velocity 
(y a = Q.llRjt) of 0.58c. The radio light curve in the 
case of a uniform density medium is relatively flat, as 
compared the wind density profile (Figure [5]), so the most 
powerful observational constraint comes from our stacked 
limit, giving a density limits in the range ncsM ^ 6 cm~ 3 



(e B = 0.1) to ricsM < 44 cm -3 (e B = 0.01). Here, our 
limits on the density scale with the assumed value for 
e B as ti-csm oc e^ ' 9 . It is noteworthy that, unlike in the 
wind case, the radio luminosity actually increases with 
time, because the decline in post-shock energy density is 
more gradual. In this uniform density case, the growth 
in the radio-emitting volume overwhelms the effects of 
gradual deceleration of the Shockwave, implying that the 
total energies in magnetic fields and relativistic electrons 
slowly grow as the SN expands. 

4. IMPLICATIONS FOR SN IA PROGENITORS 

Below, we compare these constraints on the environ- 
ment of SN 201 lfe with predictions for popular single de- 
generate scenarios, summarized in Figure [3J 
Symbiotic Systems: In the first SD scenario, the 
SN 20 life progenitor system consists of a WD in a sym- 
biotic binary system, accreting mass from a giant star 
(Figure [3]). Mass is lost from the system at a rate 
M > 10~ 8 Mr, yr" 1 with velocity y w » 30 km s" 1 
(ISeacmist fc TavlorlfT990l: iPatat et alii 2 011 at IChen et al.l 
120111 ). This model is ruled out by our EVLA observations 
and is consistent with the non -detection of a g iant star in 
pre-explosion optical imaging (|Li et a l. 2011b) and X-ray 
limits on the inverse Compton emission from t he shock- 
wave ()Horesh et al.ll2012t iMargutti et al.ll2012l ). 
White Dwarfs with Steady Nuclear Burning: 
Next, we consider a main sequence, subgiant, or helium 
star undergoing stable Roche Lobe overflow onto a WD. 
For accretion rates of M acc > 3 x 10~ 7 M w yr" 1 , the WD 
will undergo steady nuclear burning ([Shen Sz Bildstenl 
2007). We assume that a small fraction of the trans- 
ferred mass (ei oss «1%) is lost at the outer Lagrangian 
points of the system and travels outwards at at some fac- 
tor of order unity of the orbital velocity, v w ~ 100 km s _1 
(dark blue region in Figure [3]) . This mode of outer La- 
gragian mass loss with velocities up to ^600 km s _1 
has been manifested in the P-Cygni profiles of stable nu- 
clear burning WDs unde rgoing these moderate accretion 
rates (IDeufel et alJl999\ This scenario is also supported 
by the circumbina ry material observed from cataclysmic 
variable systems (Williams et al.l 120081 ) , as well as the 
steady orbital sepa rations measured i n long period cat- 
aclysmic variables (|Huang "fcYullT99fil ). Most of this pa- 
rameter space is ruled out for SN2011fe by our EVLA 
limits, in the £b = 0.1 case. 

At higher mass transfer rates, optically-thick winds in 
the vicinity of the WD are hypothesized to limit the mass 
accretion to M acc w 6 x 10~ 7 M Q yr -1 , and any ad- 
ditional mass transferred is lost from the system with 
velocity, v w ~ few x 1,000 km s _1 (light blue region 
in Figure El Hachi su"et al.1 119991 ) , consistent with out- 
flows seen in the most X- ray luminous nuclear burning 
WDs (jCowlev et al.l Il998h . This scenario is completely 
excluded by our EVLA limits if the accretion wind im- 
mediately preceded the SN explosion (see §6.1 for further 
discussion). 

Near the transition point between these moderate and 
high accretion rate regimes, spectroscopic analysis of the 
absorption profiles of luminous nuclear-burning WDs re- 
veals that the mass loss rate is ~10% the accretion rate, 
or 3 x 10 -8 Mq yr -1 , and is expell ed at a high velocity of 
v w = 3,000 km s- 1 (|Murravil2002f) . Such a high-velocity 
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Fig. 3. — Constraints on the parameter space of mass loss rate from the binary system versus the velocity of the lost material. The 
progenitor scenarios discussed in Section 4 are plotted as schematic zones in M—v w space. We indicate our 3<r limits on M /v w assuming 
eg = 0.1 (solid; corresponding to the dark blue curve in Figure 2), and the conservative case of £b = 0.01 (dashed; corresponding to the 
light red curve in Figure 2). Mass loss scenarios to the lower right of these lines would have been detected with our EVLA observations, 
and are eliminated for SN2011fe. At the lowest mass accretion rates, nova eruptions will expel the accreted material, and there will be 
insufficient mass retention toward SNe la. At high mass transfer rates > 2 X 10 -6 Mq yr — 1 , too much mass is lost though winds, so that 
the WD will not accrete sufficient mass to approach Mch- 



accretion wind would have been barely detectable by our 
radio data (assuming e_e = 0.1, see Fig. 3); nevertheless, 
the slower wind via the outer Lagragian points, albeit 
intrinsically lower in its mass loss rate, would have been 
readily detectable due to the high accretion rate in such 
a system. 

Our EVLA observations of SN2011fe therefore crit- 
ically constrain stable nuclear burning in both accre- 
tion regimes (outer Lagrangian losses and optically-thick 
winds; Figure [3]). Parallel studies based on pre-explosion 
X-ray images of the SN2011fe progenitor system only 
rule out a subset of these sys tems with the most com- 
pact and hottest photospheres (|Li et al.ll2011bt iLiu et all 
l20Tl . 

Recurrent Novae: Finally, at lower accretion rates, 
Mace w (1 - 3) x 1CT 7 M Q yr" 1 , the WD may un- 
dergo nova outbursts with a short recurrence time of 
several years. The nova shells expand with velocities 
of 1, 000 — 5, 000 km s" 1 , leaving evacuated cav ities be- 
tween the shells (|Wood-Vasev fc Sokoloskll2006[) . In the 
paradigm of synchrotron emissio n produced by th e in- 
teraction of SN shock with CSM (|Chevalierlll982b| ). ra- 



dio emission will be generated when the SN Shockwave 
crashes through a previously ejected shell. 

For SN2011fe, our EVLA observations probe a cir- 
cumbinary scale of r w 10 16 cm, thus constraining the 
presence of shells from novae with short recurrence time, 
~few years. The light green zones in Figure |4] show the 
distribution of CSM around a nova with recurrence time 
of 2 years. Models of such a system s how that 2 x 10~ 7 
M^ yr" 1 is incid ent upon the WD (jYaron et al.l 120051: 
IShen fc Bildstenl [2009) ; we assume that 1% of this trans- 
ferred material is lost to the binary to form the CSV! at 
small radii. During a nova explosion, for the adopted 
accretion rate onto a ~ 1.4 M Q WD, models imply that 
~15% of the material accreted over the last two years 
is ejected in a shell with a mass of 6 x 10~ 8 M Q (light 
green region in Figure [3]) . Adopting a shell thickness of 
~ 0.1R, the radio emission is likely short-lived (^days) 
and we estimate a ^30% probability of radio detection 
of such nova shells in our EVLA data. 

For longer recurrence times, the progenitor wind may 
enhance the local CSM density between eruptions (e.g., 
U Sco; recurrence timescale of ~10 years). The dark 
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Fig. 4. — Theoretical CSM radial density profiles are estimated for four possible SD progenitor scenarios. Red: a symbiotic progenitor 
system with M = 1 X 10~ 7 Mq yr~ 1 and v w = 30 km s _1 ; Blue: stable nuclear burning, for either outer Lagrangian losses (M = 2 X 10 — 8 
Mq yr _1 and v w = 100 km s _1 ) or optically-thick winds (M = 2 X 10 — 7 Mq yr _1 and v w = 1,000 km s —1 ); and Green: a recurrent 
nova system with short recurrence timescale (2 years; light) and longer recurrence timescale (10 years; dark). Assuming a ncsM oc r -2 
CSM density profile, we estimate the CSM densities and Shockwave radii that would correspond to our 5.9 GHz limits, for eg = 0.1 (black 
arrows) and eg = 0.01 (grey arrows). Also overplotted are limits for uniform density surroundings, again for eg = 0.1 (black dashed line) 
and eg = 0.01 (grey dashed line). 



green regions plotted in Figure [4] represent a nova with 
recurrence time of 10.3 years, modeled after U Sco. 
In the case of U Sco, a nova ejecta mass of 10 ~ 6 
Mq is e xpelled every ~10 years at a veloc i ty of 5,000 
km s " 1 (|Hachisu et al.l [2000 : ISchaeferl l20l"ot iDiaz et al.1 
120101) . Between eruptions, mass is transferred to the 
WD at a rate M^ r . = 10" 7 M yr" 1 (|Duschl et al.lll990t 
IHachisu et aL|[2000l ). A gain, we assume ei OS s ~ 1% of this 
transferred mass is lost through the outer Lagrangian 
points in a wind expelled at roughly the orbital velocity, 
v w = 100 km s _1 (dark green region in Figure[3]). The lo- 
cal density is therefore M/v w = 1 x 10~ 9 fgf w \m s- 1 ' 
and comparable to our EVLA detection limits (Fig- 
ure |U) b ut unde tectable in pre-explosion optical images 
(ILi et alJl20lTbh . 

Our assumed ejected and accreted masses for U Sco are 
intermediate in the uncertain range of estimates in the 
literature; for example, measurements of th e ejecta mass 
could be a factor of th ree higher or lower (|Evans et al.l 
2001; Dia z et ahl 120101 ) . and estimates of the accretion 
rate are model dependent. These uncertainties under- 



line our point that we can not conclusively rule out a U 
Sco-like progenitor for SN2011fe, although we can con- 
strain the large potential parameter space for U Sco-like 
progenitors. 

5. ESTIMATES OF THE BROADER 
ENVIRONMENT OF SN2011FE 

Alternatively, if SN 20 life is the result of a DD progen- 
itor system, it may have exploded into an undisturbed 
and constant density environment. As shown in Figure[TJ 
we can probe the ambient interestellar medium in the 
^0.3 kpc region surrounding the SN explosion site using 
pre-explosion H I 21 cm imaging from The H I Neary 
Galaxy Survey (THINGS: IWalter et al.ll2008l) . We mea- 
sure an integrated brightness temperature of 165.4 K km 
s _1 at the location of SN2011fe, corresponding to a H I 
column density of 3.02 x 10 20 cm~ 2 . Assuming a path 
length of 100 pc through the face-on neutral disk of M101 
and solar abundance, we roughly estimate a volume num- 
ber density of ~1 cm -3 at this location , typical of th e 
warm phase of the interstellar medium (Fcrricrc 2001). 
Our stacked EVLA limits, which imply ncsM ^ (6 — 44) 



cm 3 , are therefore consistent with expansion into the 
ambient ISM. 

Using Ha imaging as a tracer of star formation in 
M101, we can also compare the location of SN2011fc 
relative to the host galaxy's light distribution. We 
adopt the methodology for calc ulati ng a fractional flu x 
value from lFruchter et all (|2006t) and lKellv etall (|2008l) . 
We obtain an Ha n arrow-band image of the field from 
iHoopes et al.l (|2001l ) with a central wavelength of 6564 
A, and calculate the fraction of total host light in pixels 
fainter than the SN position. To determine an appro- 
priate cut-off level for the host, we create an intensity 
histogram of a 30' x 30' region centered on M101 and 
model the sky brightness distribution with a Gaussian 
profile, taking pixels with flux above 3cr to be part of 
M101. We then calculate the flux in galaxy pixels fainter 
than the location of SN2011fe and normalize by the to- 
tal galaxy flux. This gives a fractional flux value of 0.55, 
showing that SN2011fe originates from an ordinary re- 
gion in a star -forming galaxy, c onsistent with results for 
other SNe la ([Kelly et al.H2008D . This result affirms that 
SN 2011fe exploded in a typical environment for a SN la, 
and is not closely associated with a young stellar popu- 
lation (in contrast, SNe Ic and long gamma-ray bursts, 
which are thought to closely follow star formation, typi- 
cally have significantly higher fractional fluxes than SNe 
la). 

6. DISCUSSION 
6.1. Constraints on the Environments of SNe la 

Here we consider our characterization of the SN2011fc 
environment in the context of other multi-wavelength 
studies of SN la environs. 

The interaction of a SN Shockwave with surrounding 
CSM should also produce X-ray and Ha emission, along 
with the radio emission described in detail here. To date, 
searches for this emission from SNe la have yielded only 
non-detections, although they typically place weaker con- 
straints on the CSM density than radio limits (however, 
iMargutti et al.1 [20121 have recently shown that deep X- 
ray observations obtained at optical peak have the po- 
tential to surpass radio constraints). Early-time deep 
Ha observations of SN 200 lei constrain the mass loss 
rate of the progenitor to M/v w < 10 -4 ^^'g-i 
(iMattila et al.l j 2005| . For the normal SN la 2002bo, 
Hugh es et al.l (|2007l) used X-ray observations to con- 
strain M/v w < 2 x 10~ 4 1 ^ kl ^ r s -i i assuming ther- 
mal emission dominates the X-ray signal. They also re- 
yisit a previously cl aimed X-ray detection of SN 2005kc 
(jlmmler et al.1120061) . showing that this source is in fact 
not present in more carefully analyzed X-ray data. Re- 
cently, all Swift X-ray observations of SNe la to date 
have been stacked by iRussell fc Immlerl ([20121 ). yield- 
ing a non-detection of Lx < 1-7 x 10 38 erg s _1 in the 
energy range 0.2-10 keV. Radio limits on CSM around 
SN 201 lfe are consistent with deep Chandra observations 
obtained 4 days after explosion, which place an upper 

limit M /v w < 2x IP" 9 ™ °j r ~ a -L (IMargutti et al.ll20li 
see also lHoresh et al.ll2012f ) . 

Our EVLA upper limits are a strong indication that 
there is not much CSM within ~ 10 16 cm of the pro- 
genitor of SN2011fe, consistent with a uniform density 



medium of ncsM ^ 10 cm~ 3 . Similar constraints can 
be drawn from the morphology and X-ray emission of 
historical supernova remnants (SNRs) of Type la ori- 
gin. For example, X-ray studies of SNR 0509-67.5 in 
the Large Magellanic Cloud and the Tycho SNR (both 
~400 years old) imply that these SNRs are surrounded 
by uniform media of density ~0.5-l cm~ 3 (Ba denes et al.1 
120061 l200l . typical of the warm phase of the interstel- 
lar medium. Ha imaging of SN 1006 with HST shows 
that this remnant is expanding into a medium of den- 
sity 0.25-0.4 cm~ 3 , and the medium is relatively uni- 
for m, with fluctuations of ~20% on length scales of ~1 
pc (jRavmond et al.l I2007h . The youngest SNR in the 
Milky Way, Gl.9+0.3, is likely the remnant of a SN la 
that exploded ~100 years ago, and is also consistent with 
low-density surroundings (^0.03 cm~ 3 ; [ Reynold s et al.l 
l2008t IBorkowski et al.1 12010: ICarlton et al.1 120111) . This 
widespread consistency with our results is remarkable, 
as these SNRs have radii of ~ 10 19 cm, three orders of 
magnitude larger than the extent probed by our EVLA 
data. 

SNRs can also shed additional light on the 
optically-thick accre tion wind scenario, postulated by 
Hac hisu et all (|1999t) to blow from WDs that are accret- 
ing near the Eddington rate (> 6 x 10~ 7 M Q yr _1 ). We 
have already shown that our radio observations can ex- 
clude such a wind in SN 2011fe if it immediately precedes 
the SN explosion, but in other evolutionary scenarios 
there may be a delay between the production of the wind 
and the SN of - 10 5 - 10 6 years (|Han fc Podsiadlowskil 
20041). In this large cavity will surround the WD 

at the time of SN la expl osion, with a radius of ~10 pc 
and n C sM < 10~ 3 cm -3 (|Badenes et al.l l2007h . and the 
outer edge of this cavity will be marked by a dense shell 
with ncsM ~ 100 cm~ 3 . This scenario would yield radio 
non-detections at early times, as we observe, but would 
also leave a clear imprint on SNRs, wherein SNRs should 
have larger radii and faster expansion velocities than ex- 
pected for typical interstellar medium densities. Such an 
effect is not observed in a sample of seven Type la SNRs 
(jBadenes et al.|[2007D : combined with our early time ob- 
servations, this implies that accretion winds can not play 
a major role in the progenitor systems of most SNe la. 
An interaction with the warm phase of the interstellar 
medium seems to be able to explain the properties of 
most Type la SNRs. 

Another intriguing constraint on the CSM surrounding 
SNe la comes from high-resolution optical spectroscopy, 
as observations of time- variable narr ow Na I D absorp- 
tion lines in the sp e ctra of SNe la (iPatat et al.l [2 007a; 
Blondin et all 120091: iSimon et all 120091: iSternberg et al.l 
201 If ). Assuming that the variations in time are due to 
ionization by the SN light followed by gradual recom- 
bination, the distance to the absorbing material can be 
estimated at > 10 17 cm with a mass of 10~ 5 — 10~ 2 
M(T), depending on its ge ometry and dumpiness (|Chugail 
120081: ISimon et all 120091 ) . This distribution of material 
may be consistent with old nova shells from a recurrent 
nova (similar absorption features we re observed follow- 
ing the 2006 outburst of RS Oph by IPatat et ail 1201 lai 
However, U Sco did not show time- variable absorptio n 
lines during its 2010 outburst: iKafka fcW illiams 20111) . 
We note that in the case of SN 2006X (the first published 
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SN la with recognized time- variable Na I D lines), VLA 
observations were acquired and yielded a non-detection, 

constraining M/v w < 3 x 10~ 7 1 ^ ) <3 k ^ L r s -i at a distance 
of ~ 10 17 cm (Pa tat et al.ll2007al using the models de- 
scribed here). 

As we have shown, early-time radio observations can 
only rule out some of the parameter space associated 
with recurrent novae, and the observations presented 
here are sensitive to material at radii < 10 16 cm, so 
high-resolution optical spectroscopy is a complemen- 
tary strategy for probing the CSM at all scales. No 
time-variab le optical absorpti on lines are apparent in 
SN2011fe ([Patat et al.1 1201 lbl ). consistent with an ab- 
sence of CSM. We note that preliminary statistics hint 
that roughly half of SNe la occur in similarly "clean" en- 
vironments, while ~25% of SNe la in spi ral galaxies show 
narro w blue-shifted absorption lines (jSternberg et al.l 
1201 If) . If all SN la systems that show blue-shifted ab- 
sorption lines have recurrent nova progenitors, this would 
conflict with our census of recurrent novae in the local 
universe, predicting many more nov a explosions than are 
observed ijdella Valle fc LiviolH996l ). 

We emphasize, however, that the progenitors of SNe 
la may be diverse, and constraints placed on the progen- 
itor of SN2011fe will not necessarily hold for the entire 
class of SNe la. For example, while most observations of 
SN la remnants are n ot consistent with accretion winds 
([Badenes et al.ll2007D . this does not imply that all SN 
la progenitor systems do not host accretion winds. For 
one, RCW 86 may provide a counter example, as recent 
work s hows that it is likely e xpanding into a wind-blown 
cavity (I Williams et al.ll2011[ ); however, it remains uncer- 
tain if RCW 86 is a SN la or core-collapse remnant. 
Also, many SNe la do not show time-variable Na I D 
absorption lines (IPatat et al.l [2007bl: ISimon et all I2007t 
Sternberg et~aTll2011[ ) including SN2011fe (jPatat et al.l 
2011bD . This overall diversity could be due to viewing 
angle effects, differences in metallicity of the CSM, or 
a range of circumbinary environments in the progenitor 
system. 

6.2. Directly Detecting the Binary Companions of SNe 

la 

While it is possible to directly detect the massive pro- 
genitor stars of core-collapse SNe in near by galaxies us - 
ing archival pre-explosion optical images (Smartt 2009), 
both the massive WDs that produce SNe la and their 
binary comp anions are too faint to be detected in exter- 
nal galaxies ([Maoz fc MannuccH 120081 ). Indeed, no opti- 
cal source was detected pr e-explosion at the position of 
SN2011fe ([Li et al.ll2011bl ). However, in our Galaxy, we 
can search Type la SNRs for stars with unusual proper 
motions and abundances, as might be expected for the 
binary companion of the now-destroyed WD. One such 
star was proposed to be the sub-gian t star G in the Ty- 
cho SNR dRuiz-Lapuente et al.ll2004l ). but more recently 
star G's afHlilation with the SN R has been called into 
question (jKerzendorf et aLl l2009) . Additional studies are 
required, of Tycho and other SNRs in the Milky Way and 
Magellanic Clouds, to directly search for binary compan- 
ions to SNe la. 

When the SN Shockwave plows over a non-degenerate 
companion, this interaction may leave an observable 



mark upon the SN la explosion itself. If the binary 
has a relatively small separation, some of the compan- 
ion star's envelope should be entrained in the SN ejecta, 
and become d etectable in late-tim e nebular spectra as 
Ha emission (Marietta et al.l [200 0). Deep limits have 
been placed on the presence of such emission, imply- 
ing that < 0.01 M Q of hydrogen-rich material is swept 
from a companion star, and presenting a challenge to 
progenitor s cenarios where the companion is filling its 
Roche lobe ([Leonard! 12007). The interaction may also 
be detectable as an ea rly-time blue component in the 
optical SN light curve ([Kasenl 120101) ; the magnitude of 
this early component will depend on viewing angle and 
the nature of the companion, with Roche-lobe filling 
red giant stars producing the brightest signature. In 
the significant samples of SNe la now available, this 
effect has never been detected, implyin g that <20% 
of SNe la have red giant companions (|Havden et al.1 
121)101 [Bianco et al.l [20111: IBrown et all I2012D . In addi- 
tion, no blue early-time blue component was detected for 
SN2011fe, placing strong constraints on possible Roche- 
lobe filling co mpanions (<1 M B red gian t, or <3 M^ 
main sequence: IBrown et al.ll201 it : see also lNugent et al.l 
[20Tl1 : IBloom et al.N2012l )~ 

7. CONCLUSIONS 

Using deep EVLA radio observations, we have shown 
that the progenitor system of SN2011fc did not host 
a red giant secondary or a wind-producing high accre- 
tion rate. We also rule out a significant fraction of 
the parameter space for stably burning WD or recur- 
rent nova progenitors — constraints that are bolstered by 
pre-explosion Chandra and HST imaging and a lack of 
an early-time blue component in the optical light curve. 

Our study may be viewed as consistent with the grow- 
ing body of evidence favoring DD progenitors, which 
includes m easurements of th e delay time distribution 
of SNe la ([Maoz et al.l 120101 ) and comp utational suc- 
cesse s with theoretical DD explosions (jPakmor et al.1 
l2010f) . However, a recent preliminary model of a WD- 
WD mer ger suggest s the p resence of dense circumbinary 
medium IShen et ail ([201 1|) . which could be constrained 
by these and future radio observations. A thorough treat- 
ment of this scenario would require a re-evaluation of the 
SN Shockwave dynamics, as it is predicted that a dense 
envelope of material will surround the primary WD at 
small radius (^0.1 M Q at < 10 13 cm), which would slow 
the Shockwave. At larger radius, the WD is surrounded 
by a pcsm oc r~ 2 wind profile, expelling roughly ~0.1 
Mq over 10 4 years at a velocity of ^100 km s , produc- 
ing a relatively high M/v w sa 10~ 5 s -i ■ It is the 
interaction between this wind and the SN Shockwave that 
would likely produce detectable radio emission, but more 
detailed models of the density profile and mass loss in DD 
systems are needed to understand the velocity evolution 
of the Shockwave and the predicted synchrotron signal. 

More exotic scenarios, like spin- up/spi n-down mod- 
els (|Justhaml 120111 : iDi Stefano et al.ll201lD or the core- 
degenerate scenario ([Ilkov fc Soken 120121 ) . with a signifi- 
cant delay between mass loss and explosion (> 10 5 years) 
are also not excluded. While our EVLA upper limits 
place the most stringent constraints to date on the CSM 
density around a SN la, and our analysis indicates that 
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most SD progenitor models are ruled out, detailed stud- 
ies (theoretical and observational) of the circumbinary 
environments surrounding accreting WDs are required 
to further shed light on their connection to Type la SNe. 
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